A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.